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We present a theoretical analysis of the connection between classical polarization optics and quan- 
tum mechanics of two-level systems. First, we review the matrix formalism of classical polarization 
optics from a quantum information perspective. In this manner the passage from the Stokes- Jones- 
Mueller description of classical optical processes to the representation of one- and two-qubit quantum 
operations, becomes straightforward. Second, as a practical application of our classical-ns-quantum 
formalism, we show how two-qubit maximally entangled mixed states (MEMS), can be generated 
by using polarization and spatial modes of photons generated via spontaneous parametric down 
conversion. 



PACS numbers: 03.65.Ud, 03.67.Mn, 42.25.Ja 



I. INTRODUCTION 

Quantum computation and quantum information have 
been amongst the most popular branches of physics in 
the last decade One of the reasons of this success 
is that the smallest unit of quantum information, the 
qubit, could be reliably encoded in photons that are easy 
to manipulate and virtually free from decohcrcnce at op- 
tical frequencies @, Q. Thus, recently, there has been a 
growing interest in quantum information processing with 
linear optics [1, [U, [g, 0] and several techniques to gener- 
ate and manipulate optical qubits have been developed 
for different purposes ranging from, e.g., teleportation 
@> @, to quantum cryptography to quantum mea- 
surements of qubits states [l(| and processes [ll|, etc. 
In particular, Kwiat and coworkers [H, [l3| were able 
to create and characterize arbitrary one- and two-qubit 
states, using polarization and frequency modes of pho- 
tons generated via spontaneous parametric down conver- 
sion (SPDC) [3. 

Manipulation of optical qubits is performed by means 
of linear optical instruments such as half- and quarter- 
wave plates, beam splitters, polarizers, mirrors, etc., and 
networks of these elements. Each of these devices can 
be thought as an object where incoming modes of the 
electromagnetic fields are turned into outgoing modes by 
a linear transformation. From a quantum information 
perspective, this transforms the state of qubits encoded 
in some degrees of freedom of the incoming photons, ac- 
cording to a completely positive map £ describing the 
action of the device. Thus, an optical instrument may 
be put in correspondence with a quantum map and vice 
versa. Such correspondence has been largely exploited 
0, 0, El Ell and stressed [H, E3 by several authors. 
Moreover, classical physics of linear optical devices is a 
textbook matter [18|, and quantum physics of ele- 
mentary optical instruments has been studied extensively 
[20| . as well. However, surprisingly enough, a systematic 
exposition of the connection between classical linear op- 
tics and quantum maps is still lacking. 

In this paper we aim to fill this gap by presenting a de- 
tailed theory of linear optical instruments from a quan- 



tum information point of view. Specifically, we estab- 
lish a rigorous basis of the connection between quantum 
maps describing one- and two-qubit physical processes 
operated by polarization-affecting optical instruments, 
and the classical matrix formalism of polarization optics. 
Moreover, we will use this connection to interpret some 
recent experiments in our group [21| . 

We begin in Section II by reviewing the classical the- 
ory of polarization-affecting linear optical devices. Then, 
in Section III we show how to pass, in a natural manner, 
from classical polarization-affecting optical operations to 
one-qubit quantum processes. Such passage is extended 
to two-qubit quantum maps in Section IV. In Section V 
we furnish two explicit applications of our classical- vs- 
quantum formalism that illustrate its utility. Finally, in 
Section V we summarize our results and draw the con- 
clusions. 



II. CLASSICAL POLARIZATION OPTICS 

In this Section we focus our attention on the descrip- 
tion of non-image-forming polarization-affecting optical 
devices. First, we shortly review the mathematical for- 
malism of classical polarization optics and establish a 
proper notation. Second, we introduce the concepts of 
Jones and Mueller matrices as classical maps. 



A. Polarization states of light beams 

Many textbooks on classical optics introduce the con- 
cept of polarized and unpolarized light with the help of 
the Jones and Stokes-Mueller calculi, respectively [18l |. 
In these calculi, the description of classical polarization 
of light is formally identical to the quantum description 
of pure and mixed states of two-level systems, respec- 
tively [221 ]. In the Jones calculus, the electric field of a 
quasi- monochromatic polarized beam of light which prop- 
agates close the ^-direction, is represented by a complex- 
valued two-dimensional vector, the so-called Jones vec- 
tor E e C 2 : E = E x + Eiy, where the three real- 
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valued unit vectors {x, y, z} define an orthogonal Carte- 
sian frame. The same amount of information about the 
state of the field is also contained in the 2x2 matrix J 
of components Jy = EiE* , (i,j = 0, f), which is known 
as the coherency matrix of the beam [19[ . The matrix J 
is Hcrmitcan and positive scmidcfinitc 

jt = J (v,Jv) = |(v ! E)| 2 >0, (1) 

= Yli=o u t v i denotes the ordi- 
Further, J has the projection 



where v £ C , and (u,v) 
nary scalar product in C 2 
property 



J = JTrJ, 



(2) 



and its trace equals the total intensity of the beam: 
TrJ = |_E | 2 + \E\\ 2 . If we choose the electric field units 
in such a way that TrJ = 1, then J has the same proper- 
ties of a density matrix representing a two-level quantum 
system in a pure state. In classical polarization optics the 
coherency matrix description of a light beam has the ad- 
vantage, with respect to the Jones vector representation, 
of generalizing to the concept of partially polarized light. 
Formally, the coherency matrix of a partially polarized 
beam of light is characterized by the properties |Tj), while 
the projection property ([2]) is lost. In this case J has the 
same properties of a density matrix representing a two- 
level quantum system in a mixed state. Coherency matri- 
ces of partially polarized beams of light may be obtained 
by tacking linear combinations ^2 n wnJn of coherency 
matrices Jn of polarized beams (all parallel to the same 
direction z), where the index N runs over an ensemble 
of field configurations and wn > 0. The degree of polar- 
ization (DOP, denoted P) of a partially polarized beam 
is defined by the relation 

Det J = (TrJ) 2 (1 - P 2 )/4. (3) 

For a polarized beam of light, projection property © 
implies Det J = and P = 1, otherwise < P < 1. 
It should be noted that the off-diagonal elements of the 
coherency matrix are complex-valued and, therefore, not 
directly observables. However, as any 2x2 matrix, J can 
be written either in the Pauli basis X a : 



Xn = 
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(4) 



or in the Standard basis Y a : 
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(5) 



as 



^ 3 3 

j = ^ Y XaXa = Y yp Y p> 



(6) 



where x a = Tr(X a J) elj ? = Tr(Y^J) £ C and, from 
now on, all Greek indices a, /3, fi, u, . . . , take the values 
0, 1,2,3. The four real coefficients x a , called the Stokes 
parameters [23j of the beam, can be actually measured 
thus relating J with observables of the optical field. For 
example, xq = TrJ represents the total intensity of the 
beam. Conversely, the four complex coefficients yp are 
not directly measurable but have the advantage to fur- 
nish a particularly simple representation of the matrix J 
since y Q = J 00 , yi = Jn, V% = Jio, 2/3 = Jn- The two 
different representations x a and yp are related via the 
matrix 
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(7) 



such that x a = Vapyp, where V a p = Ti:{X a Yi3), and 
V^V = 2/4 = VV\ where I4 is the 4x4 identity matrix. 



B. Polarization-transforming linear optical 
elements 



When a beam of light passes through an optical system 
its state of polarization may change. Within the con- 
text of polarization optics, a polarization- affecting linear 
optical istrument is any device that performs a linear 
transformation upon the electric field components of an 
incoming light beam without affecting the spatial modes 
of the field. Half- and quarter- wave plates, phase shifters, 
polarizers, are all examples of such devices. The class of 
polarization-affecting linear optical elements comprises 
both non-depolarizing and depolarizing devices. Roughly 
speaking, a non- depolarizing linear optical element trans- 
forms a polarized input beam into a polarized output 
beam. On the contrary, a depolarizing linear optical ele- 
ment transforms a polarized input beam into a partially 
polarized output beam [241 ]. A non-depolarizing device 
may be represented by a classical map via a single 2x2 
complex- valued matrix T, the Jones matrix [1 81 ] . such 
that 



TE, 



(8) 



for polarized input beams or, for light beams with arbi- 
trary degree of polarization: 



(9) 



(3=0 



In this paper we consider only passive (namely, non- 
amplifying) optical devices for which the relation 
TrJ out < TrJj n holds. There exist two fundamental 
kinds of non-depolarizing optical elements, namely re- 
tarders and diattenuators; any other non-depolarizing 
element can be modelled as a retarder followed by a 
diattenuator (25j. A retarder (also known as birefrin- 
gent element) changes the phases of the two components 
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of the electric-field vector of a beam, and may be rep- 
resented by a unitary Jones matrix Ty. A diattenu- 
ator (also known as dichroic element) instead changes 
the amplitudes of components of the electric-field vector 
(polarization-dependent losses), and may be represented 
by a Hermitean Jones matrix Tjj. 

Let 7nd denotes a generic non-depolarizing device rep- 
resented by the Jones matrix T, such that J ul — > J ou t = 
TJi n T^. We can rewrite explicitly this relation in terms 
of components as 



\Jo\it)ij — TikTji(Ji n )kl, 



(10) 



where, from now on, summation over repeated indices 
is understood and all Latin indices i, j, k, l,m,n, . . . take 
the values and 1. Since T lk T* = {T®T*) tjM = M ljM 
we can rewrite Eq. (TIT)]) as 

(11) 



{Jontjij = M.ij,kl(Jin)kl, 



where M. = T ®T* is a 4 x 4 complex-valued matrix 
representing the device 7nd, and the symbol ® denotes 
the ordinary Kronecker matrix product. M is also known 
as the Mueller matrix in the Standard matrix basis [26[ 
and it is simply related to the more commonly used real- 
valued Mueller matrix M [18J via the change of basis 
matrix V: 



M= -VMVl 



(12) 



For the present case of a non-depolarizing device, M is 
named as Mueller- Jones matrix. From Eqs. ((6j ITTj) it 
readily follows that we can indifferently represent the 
transformation operated by 7nd either in the Standard 
or in the Pauli basis as 
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0=0 
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,9=0 



M a/ 3x%, (13) 



respectively. 

With respect to the Jones matrix T, the Mueller matri- 
ces M. and M have the advantage of generalizing to the 
representation of depolarizing optical elements. Mueller 
matrices of depolarizing devices may be obtained by tak- 
ing linear combinations of Mueller-Jones matrices of non- 
depolarizing elements as 



M = J^PaMa = ^PaTa ® T* A , 



(14) 



where pa > 0. Index A runs over an ensemble (either de- 
terministic 27] or stochastic 28]) of Mueller- Jones matri- 
ces Ma — Ta®T\, each representing a non-depolarizing 
device. The real-valued matrix M corresponding to M. 
written in Eq. (|T3|), can be easily calculated by using Eq. 
([12)1 that it is still valid 26] . In the current literature M 
is often written as \2M 



M = 



Moo d T 
p W 



(15) 



where pg I 3 , d € R 3 , are known as the polarizance vec- 
tor and the diattenuation vector (superscript T indicates 
transposition) , respectively, and W is a 3 x 3 real- valued 
matrix. Note that p is zero for pure depolarizers and pure 
retarders, while d is nonzero only for dichroic optical ele- 
ments |25l j . Moreover, W reduces to a three-dimensional 
orthogonal rotation for pure retarders. It the next Sec- 
tion, we shall show that if we choose Moo = 1 (this can 
be always done since it amounts to a trivial polarization- 
independent renormalization) , the Mueller matrix of a 
non-dichroic optical element (d = 0), is formally identi- 
cal to a non-unital, trace-preserving, one-qubit quantum 
map (also called channel) [29j]. If also p = (pure de- 
polarizers and pure retarders), then M is identical to a 
unital one-qubit channel (as defined, e.g., in [l[). 



III. FROM CLASSICAL TO QUANTUM MAPS: 
THE SPECTRAL DECOMPOSITION 



An important theorem in classical polarization optics 
states that any linear optical element (either determinis- 
tic or stochastic) is equivalent to a composite device made 
of at most four non-depolarizing elements in parallel [30l ] . 
This theorem follows from the spectral decomposition of 
the Hermitean positive semidefinite matrix H [311 ] as- 
sociated to Ai. In this Section we shortly review such 
theorem and illustrate its equivalence with the Kraus de- 
composition theorem of one-qubit quantum maps [l|. 

Given a Mueller matrix A4, it is possible to built a 
4x4 Hermitean positive semidefinite matrix H — H(A4) 
by simply reshuffling [32| the indices of A4 : 



Hij,ki = Mikji = ^pa(Ta)ij(T^ 



AJkl, 



(16) 



where the last equality follows from Eq. (fLT| . Equiv- 
alently, after introducing the composite indices a — 
2i + j, [3 — 2k + I, we can rewrite Eq. (fTI)]) as H a p = 
*}2,APA{TA)a{TX)f3- In view of the claimed connection 
between classical polarization optics and one-qubit quan- 
tum mechanics, it worth noting that H is formally identi- 
cal to the dynamical (or Choi) matrix, describing a one- 
qubit quantum process [33]. The spectral theorem for 
Hermitean matrices provides a canonical (or spectral) 
decomposition for H of the form [34[ 



H, 



a/3 



M=0 



A (U (u /1 ) Q (u*) /3 , 



(17) 



where A M > are the non-negative eigenvalues of H, 
and {u M } = {uo, ui, 112, 113} is the orthonormal basis 
of eigenvectors of H: Hu^ = A^u^. Moreover, from a 
straightforward calculation it follows that: X)^=o 'V = 
2Mqo [26j . If we rearrange the four components of each 
eigenvector u M to form a 2 x 2 matrices T M defined as 



(u M ) (u p )i 
(u M ) 2 KJ3 



(18) 
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we can rewrite Eq. Q2|) as H af} = £^ A Al (T tI ) ct (T*) ( g. 
Since Eq. (fl~8f can be rewritten as (T^)y = (u ft ) et =2i+j, 
we can go back from Greek to Latin indices and rewrite 
Eq. dTTJ) as 

3 3 

n ijtkl = *Ma(Tpki = E v(^® t ;w- (is) 

|U=0 M=0 

Finally, from the relation above and using Eq. (|16j) . we 
obtain 

3 

w = E a ^® t ;- ( 2o ) 

Equation (|20|1 represents the content of the decomposi- 
tion theorem in classical polarization optics, as given by 
Cloude [3(| HH . It implies, via Eq. (fTTj) , that the most 
general operation that a linear optical device can perform 
upon a beam of light can be written as 

3 

J'm - * ^out = A^T^JinT^, (21) 

where the four Jones matrices represent four different 
non-depolarizing optical elements. 

Since A M > 0, Eq. pTjl is formally identical to the 
Kraus form Q of a completely positive one-qubit quan- 
tum map £. Therefore, because of the isomorphism be- 
tween J and p [22| . when a single photon encoding a 
polarization qubit (represented by the 2x2 density ma- 
trix /0i n ), passes through an optical device classically de- 
scribed by the Mueller matrix M = A M T M ® T*, its 
state will be transformed according to 

3 

Pin — ► Pout oc ^A^T^Tt, (22) 

where the proportionality symbol "oc" accounts for a pos- 
sible renormalization to ensure Trp out = 1. Such renor- 
malization is not necessary in the corresponding classical 
equation (|21[) since Tr J out is equal to the total intensity 
of the output light beam that does not need to be con- 
served. Note that by using the definition (|2"0|) we can 
rewrite explicitly Eq. (f2"2"| as 

Pout,ij OC p ou t,ij = ■Mij,klPin,kh (23) 

where (p)ij = (i\p\j) are density matrix elements in the 
single-qubit standard basis i G {0,1}, and p out is 

the un-normalized single-qubit density matrix such that 
Pout = Pout/Trp out . From Eqs. (HI5J> and Eq. J23J), it 
readily follows 

Trp out = M 00 + Moi(pi„ i0 i + Pin.io) 

+iM 2(pi,i,01 - Pin,io) 

+M 03 (pi„,oo -Pin.il), (24) 



where we have assumed Trpj n = 1. The equation above 
shows that M represents a trace-preserving map only if 
Moo = 1 and d T = (M i, M 02 , M 03 ) = (0,0,0), namely, 
only if M. describes the action of a non-dichroic optical 
instrument. In addition, if p ln represents a completely 
mixed state, that is if p m = Xq/2, then from Eq. (|23p it 
follows: 

1 3 

Pout = -j^P^X^, (25) 

were we have defined po = Moo and (pi,P2,P3) = P is the 
polarizance vector. Equation (|25[) shows that in this case 
Trp ou t = M 00 , and p out = Pout/Mao ^ X /2 if p ^ 0, 
that is, the map represented by M. (or, M) is unital only 
if p = 0. 

By writing Eqs. l|2U25j) we have thus completed the 
review of the analogies between linear optics and one- 
qubit quantum maps. In the next Section we shall study 
the connection between classical polarization optics and 
two- qubit quantum maps. 

IV. POLARIZATION OPTICS AND 
TWO-QUBIT QUANTUM MAPS 

Let us consider a typical SPDC setup where pairs of 
photons are created in the quantum state p along two 
well defined spatial modes (say, path A and path B) of 
the electromagnetic field, as shown in Fig. 1. Each pho- 




FIG. 1: (Color online) Layout of a typical SPDC experimental 
setup. An optically pumped nonlinear crystal, emits photon 
pairs that propagate along path A and B through the scatter- 
ing devices Ta and Tb, respectively. Scattered photons are 
detected in coincidence by detectors Da and Db that permit 
a tomographically complete two-photon polarization state re- 
construction. 

ton of the pair encodes a polarization qubit and p can 
be represented by a 4 x 4 Hermitean matrix. Let Ta 
and Tb be two distinct optical devices put across path 
A and path B, respectively. Their action upon the two- 
qubit state p can be described by a bi-local quantum 
map p — > Ea ® £b[p\ [Hi- A sub-class of bi-local quan- 
tum maps occurs when either Ta or Tb is not present 
in the setup, then either Ea — T or Eb = T, respec- 
tively, and the corresponding map is said to be local. In 
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the above expressions X represents the identity map: It 
does not change any input state. When a map is local, 
that is when it acts on a single qubit, it is subjected to 
some restrictions. This can be easily understood in the 
following way: For definiteness, let assume £b = T so 
that the local map £ can be written as £[p] — £a ®I[p\- 
Let Alice and Bob be two spatially separated observer 
who can detect qubits in modes A and B, respectively, 
and let p and ps denote the two-qubit quantum state be- 
fore and after Ta, respectively. In absence of any causal 
connection between photons in path A with photons in 
path B, special relativity demands that Bob cannot de- 
tect via any type of local measurement the presence of 
the device Ta located in path A. Since the state of each 
qubit received by Bob is represented by the reduced den- 
sity matrix p B = Tr| a{ps)i the locality constraint can be 
written as 

pf = p B . (26) 

We can write explicitly the map £a ® T as a Kraus 
operator-sum decomposition [l[ 

3 

p H-> p £ OC 

where, from now on, the symbol / denotes the 2x2 
identity matrix and {A^} is a set of four 2x2 Jones 
matrices describing the action of Ta- Then, Eq. (|26]) 
becomes 

3 

k,l fi=0 k 

which implies the trace-preserving condition on the local 
map £a <8> T: 

3 

X^A+A^oc/. (29) 

fi=0 

Local maps that do not satisfy Eq. (12"9"|) are classified as 
non-physical. In this Section we show how to associate 
a general two-qubit quantum map £[p] = £ a <8> £b[p] to 
the classical Mueller matrices A4 A and M B describing 
the optical devices Ta and Tb , respectively. Surprisingly, 
we shall find that do exist physical linear optical devices 
(dichroic elements) that may generate non-physical two- 
qubit quantum maps (37j . 

Let denotes with \ij) = \i) ® \j), i,j S {0, 1} the two- 
qubit standard basis. A pair of qubits is initially prepared 
in the generic state p = p ij:k i\ij)(kl\ = pf ktjl \i)(k\ <g) 
where superscript R indicates reshuffling of the indices, 
the same operation we used to pass from M to H: 
pf k - t = pij t ki = (ij\p\kl). P is transformed under the 
action of the bi-local linear map £[p] = £a® £b[p] into 
the state 

p £ = £ A ®£ B \p) o:^\ l i\ v {A il ®B v )p(Al®Bt), (30) 



where {A^} and {B u } are two sets of 2 x 2 Jones matrices 
describing the action of Ta and Tb, respectively. From 
Eq. (|3"0|) we can calculate explicitly the matrix elements 
(ij\p£\kl) = (ps)ij^ki in the two-qubit standard basis: 

(pe)ij,ki oc A^(A (li ) im (A*)/ S pp* png X v (B v ) jn (B*)i q 

= M A \A B o R 

J ik,mp J l jl,nqrmp,nq' 

(31) 

where summation over repeated Latin and Greek indices 
is understood. Since by definition (pe)ij^ki = (Ps)ik,ji we 
can rewrite Eq. (|3ip using only Greek indices as 

{p*U « M%M% P % = (M A ® M B ) ap ^ p%, 

(32) 

where summation over repeated Greek indices is again 
understood. Equation (j3"2"| relates classical quantities 
(the two Mueller matrices M A and M. B ) with quan- 
tum ones (the input and output density matrices p R and 
p R , respectively). Moreover, it is easy to see that Eq. 
(|52")) is the two-qubit quantum analogue of Eq. (|13p . 
In fact, if we introduce the 16 x 16 two-qubit Mueller 
matrix M. = M A ® A4 B , and the input and output 
two-qubit Stokes parameters in the standard basis de- 
fined as: yf =itl+v = (p R )^, y™4a+/3 = (Ps) a p, where 
a, b G {0, . . . , 15}, then we can write Eq. ([32)) as 

15 

y r«^AW n , (33) 

6=0 

which is formally identical to Eq. (H]). Thus, Eq. Q 
realizes the connection between classical polarization op- 
tics and two-qubit quantum maps. 

An important case occurs when £b = T =>• M. B = I4 
and Eq. (f32|) reduces to 

pf oc M A p R . (34) 

Equation (|34[) illustrates once more the simple relation 
existing between the classical Mueller matrix M A and 
the quantum state ps- 

With a typical SPDC setup it is not difficult to pre- 
pare pairs of entangled photons in the singlet polariza- 
tion state. Via a direct calculation, it is simple to show 
that when p represents two qubits in the singlet state 
p s = \{X Q ®X Q -X l ®X 1 -X 2 ®X 2 -X ? ,® X 3 ) and 
M A is normalized in such a way that M A = 1, then the 
proportionably symbol in the last equation above can be 
substituted with the equality symbol: 

P§ = Mp R => p £ = (Mp R ) R , (35) 

where, from now on, we write M for M A to simplify the 
notation. Note that this pleasant property is true not 
only or the singlet but for all four Bell states [l| , as well. 
Equation (|35p has several remarkable consequences: Let 
M denotes the real-valued Mueller matrix associated to 
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1. Then, the following results 



M and assume Moo 
hold: 

Tr(pf) = Tr(MM T )/4, (36) 
Tt\ A (fie) = (A + D)+M i{B + C) 

+zM 02 (B-C) + M 03 (A - D) , (37) 

where p,*; = (Mp R ) R is the un- normalized output density 
matrix. Equation (|3T[) is more general than Eq. (|36p . 
since it holds for any input density matrix p and not 
only for the singlet one p s . In addition, in Eq. (I37p we 
wrote the input density matrix p in a block-matrix form 
as 

"A fl 
C £> 



P = 



(38) 



where A, B, C = B\ and D are 2x2 sub-matrices 
and A + D = Tr\ a(p)- Equation p6|) shows that the 
degree of mixedness of the quantum state pg is in a 
one-to-one correspondence with the classical depolariz- 
ing power [24| of the device represented by M. Fi- 
nally, Eq. (|3ll), together with Eqs. I|15I26[) . tells us 
that the two-qubit quantum map Eq. (|35[) is trace- 
preserving only if the device is not dichroic, namely only 
if d T = (M 01 ,M 02 ,M 03 ) = (0,0,0). This last result 
shows that despite of their physical nature (think of, e.g., 
a polarizer), dichroic optical elements must be handled 
with care when used to build two-qubit quantum maps. 
We shall discuss further this point in the next Section. 

Before concluding this Section, we want to point out 
the analogy between the 16 x 16 Mueller matrix M = 
M A <g> M B associated to a bi- local two-qubit quantum 
map, and the 4x4 Mueller- Jones matrix M — T®T* rep- 
resenting a non-depolarizing device in a one-qubit quan- 
tum map. In both cases the Mueller matrix is said to 
be separable. Then, in Eq. (|14p we learned how to build 
non-separable Mueller matrices representing depolarizing 
optical elements. By analogy, we can now build non- 
separable two-qubit Mueller matrices representing non- 
local quantum maps, as 

M = ^w AB M A ® M b , (39) 

A,B 

where wab > 0, wab ^ wa x wb, and indices A,B run 
over two ensembles of arbitrary Mueller matrices M 
and M B representing optical devices located in path A 
and path B, respectively. 



V. APPLICATIONS 

In this Section we exploit our formalism, by apply- 
ing it to two different cases. As a first application, we 
build a simple phenomcnological model capable to ex- 
plain certain of our recent experimental results (2lj about 
scattering of entangled photons. The second application 
consists in the explicit construction of a bi-local quan- 
tum map generating two-qubit MEMS states. A realistic 
physical implementation of such map is also given. 



A. Example 1: A simple phenomenological model 

In Ref. [2l[, by using a setup similar to the one 
shown in Fig. 1, we have experimentally generated en- 
tangled two-qubit mixed states that lie upon and below 
the Werner curve in the linear entropy-tangle plane [38| . 
In particular, we have found that: (a) Birefringcnt scat- 
terers always produce generalized Werner states of the 
form pew = V ® I pwV^ ® I , where pw denotes ordinary 
Werner states (39j , and V represents an arbitrary unitary 
operation; (b) Dichroic scatterers generate sub- Werner 
states, that is states that lie below the Werner curve in 
the linear entropy-tangle plane. In both cases, the input 
photon pairs were experimentally prepared in the polar- 
ization singlet state p s . In this subsection we build, with 
the aid of Eq. (|35p . a phenomenological model explaining 
both results (a) and (b). 

To this end let us consider the experimental setup rep- 
resented in Fig. 1. According to the actual scheme used 
in Ref. [2l[, where a single scattering device was present, 
in this Subsection we assume Tg = X, so that the re- 
sulting quantum map is local. The scattering element 
Ta inserted across path A can be classically described by 
some Mueller matrix M. In Ref. (25j, Lu and Chipman 
have shown that any given Mueller matrix M can be 
decomposed in the product 



M = MdMbMa, 



(40) 



where A^a, Mb, and Md are complex- valued Mueller 
matrices representing a pure depolarizer, a retarder, and 
a diattenuator, respectively. Such decomposition is not 
unique, for example, M = MaMdMb is another valid 
decomposition [40| . Of course, the actual values of Ma, 
Mb, and Md depend on the specific order one chooses. 
However, in any case they have the general forms given 
below: 



M A = 



2 " 

n a + t 

r\ a—b a+b 
2 ° 

^ 





-b 

2 , 
a+b 

2 





1-c 
2 




1+c 
2 



Mb 
M d 



= T H ®T, 



H- 



(41) 

(42) 
(43) 



where a,b,c G M, and Tjj, Th are the unitary and 
Hermitean Jones matrices representing a retarder and 
a diattenuator, respectively. Actually, the expression 
of A^ a g iven in Eq. (|41| is not the most general pos- 
sible [25(, but it is the correct one for the representa- 
tion of pure depolarizers with zero polarizance, such as 
the ones used in Ref. [13] • Note that although M b 
and Md are Mueller- Jones matrices, Ma is not. When 
a = b = c = p : p£ [0,1] the depolarizer is said to be 
isotropic (or, better, polarization-isotropic). This case 
is particularly relevant when birefringence and dichroism 
are absent. In this case Mb = h = Md, and Eq. (|40|) 
gives M = Ma- Thus, by using Eq. (|4"Tj) we can calcu- 
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late Ma(p) and use it in Eq. ([35]) to obtain 



Ps = PPs H -r—h 



Pw, 



(44) 



that is, we have just obtained a Werner state: = p\\-\ 
Thus, we have found that a local polarization-isotropic 
scatterer acting upon the two-qubit singlet state, gener- 
ates Werner states. 

Next, let us consider the cases of birefringent (re- 
tarders) and dichroic (diattenuators) scattering devices 
that we used in our experiments. In these cases the total 
Mueller matrices M. of the devices under consideration, 
can be written as A4 = Mz-Ma, where either Z = B 
or Z = D, and A4a = Ma{p) represents a polarization- 
isotropic depolarizer. For definiteness, let consider in de- 
tail only the case of a birefringent scatterer, since the 
case of a dichroic one can be treated in the same way. In 
this case 



M b Ma(p) = J2 X v(p) t uTk 



(45) 



and, as result of a straightforward calculation, Ao = (1 + 
3p)/2, Ai = A 2 = A 3 = (l-p)/2, T p = while T v 

is an arbitrary unitary 2x2 Jones matrix representing a 
generic retarder. For the sake of clarity, instead of using 
directly Eq. ([55)1 . we prefer to rewrite Eq. ([50)1 adapted 
to this case as 



we generated 10 4 states pg from Eq. ([47]) . by randomly 
generating (with uniform distributions) the four param- 
eters p, do, d\, and 9 in the ranges: p,do,d\ € [0,1], 
9 G (0, 27r]. The numerical simulation shows that a local 




0.2 0.4 0.6 0.! 

Linear entropy, Sl 

FIG. 2: Numerical simulation from our phenomenological 
model qualitatively reproducing the behavior of a dichroic 
scattering system. The gray region represents unphysical 
states and it is bounded from below by MEMS (dashed curve). 
The lower continuous thick curve represents Werner states. 

dichroic scatterer may generate sub- Werner two-qubit 
states, that is states located below the Werner curve in 
the tangle-linear entropy plane. The qualitative agree- 
ment between the result of this simulation and the exper- 
imental findings shown in Fig. 3 of Ref. 21| is evident. 



P£ = E X » (P) ® I) Ps (Tfil ® 7) 



T V ®I 



J2K(P){T^I) Ps (Tt®j) 

_M=0 



Tjj ® I 



(46) 



where Eq. (|44p has been used. Equation Eq. ([4^)1 clearly 
shows that the effect of a birefringent scatterer is to gen- 
erate what we called generalized Werner states, in full 
agreement with our experimental results [2l| . 

The analysis for the case of a dichroic scatterer can be 
done in the same manner leading to the result 



p £ <xp £ = T H ® IpwTjj <g> I, 



(47) 



where Tjj is a 2 x 2 Hcrmitcan matrix representing a 
generic diattenuator (l8j : 



T, 



H 



do cos 9 2 + di sin 6 
(do — di ) cos 9 sin 



2 (do — di) cos 6> sin ( 
9 di cos 9 2 + d sin 6 



,(48) 



Discussion 

It should be noticed that while we used the equality 
symbol in writing Eq. (|46p , we had to use the proportion- 
ality symbol in writing Eq. (147|) . This is a consequence of 
the Hermitean character of the Jones matrix Tjj that gen- 
erates a non-trace-preserving map. In fact, in this case 
from M = M D M A (p), where M D = [VT H <g> T^)/2 
and M A (p) = [VM A (p)V^}/2 [see Eq. pj])]. we obtain 
Tr(p f ) = (dl + d\)/2 f 1. Moreover, Eq. (37]) gives 



TrU(p£ 



2 



d§ 



d\ \ Xisin26»-r-X 3 cos26' 



dg 



d? 



(49) 



where pg — pg /Tr(pg). This result is in contradiction, 
for do d\, with the locality constraint expressed by Eq. 
PB| which requires 



(50) 



where di € [0,1], are the diattenuation factors, while 
9 E (0, 27r] gives the direction of the transmission axis 
of the linear polarizer to which Th reduces when either 
do = or di = 0. Figure 2 reports, in the tangle- linear 
entropy plane, the results of a numerical simulation were 



As we already discussed in the previous Section, only 
the latter result seems to be physically meaningful since 
photons in path B, described by pf , cannot carry infor- 
mation about device Ta which is located across path A. 
On the contrary, Eq. ([4"9")l shows that p^ is expressed in 
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terms of the four physical parameters p, do, d± and 9 that 
characterize T4 . Is there a contradiction here? 

In fact, there is none! One should keep in mind that 
Eq. (|49p expresses the one-qubit reduced density matrix 
p B that is extracted from the two-qubit density matrix 
ps after the latter has been reconstructed by the two 
observers Alice and Bob by means of nonlocal coinci- 
dence measurements. Such matrix contains information 
about both qubits and, therefore, contains also informa- 
tion about T4. Conversely, p B = X /2 in Eq. (fST)|) . is 
the reduced density matrix that could be reconstructed 
by Bob alone via local measurements before he and Alice 
had compared their own experimental results and had 
selected from the raw data the coincidence counts. 

From a physical point of view, the discrepancy be- 
tween Eq. (|4"9")> and Eq. ([50)1 is due to the polarization- 
dependent losses (that is, do d\) that characterize 
dichroic optical devices and it is unavoidable when such 
elements are present in an experimental setup. Actually, 
it has been already noticed that a dichroic optical element 
necessarily performs a kind of post-selective measure- 
ment [16j . In our case coincidence measurements post- 
select only those photons that have not been absorbed 
by the dichroic elements present in the setup. However, 
since in any SPDC setup even the initial singlet state 
is actually a post-selected state (in order to cut off the 
otherwise overwhelming vacuum contribution) , the prac- 
tical use of dichroic devices does not represent a severe 
limitation for such setups. 



B. Example 2: Generation of two-qubit MEMS 
states 



In the previous subsection we have shown that it is pos- 
sible to generate two-qubit states represented by points 
upon and below the Werner curve in the tangle-linear 
entropy plane, by operating on a single qubit (local op- 
erations) belonging to a pair initially pre par ed in the 
entangled singlet state. In another paper [37| we have 
shown that it is also possible to generate MEMS states 
(see, e.g., [H, Hl[ and references therein), via local op- 
erations. However, the price to pay in that case was 
the necessity to use a dichroic device that could not be 
represented by a "physical" , namely a trace-preserving, 
quantum map. In the present subsection, as an example 
illustrating the usefulness of our conceptual scheme, we 
show that by allowing bi-local operations performed by 
two separate optical devices Ta and Tb located as in Fig. 
1, it is possible to achieve MEMS states without using 
dichroic devices. 

To this end, let us start by rewriting explicitly Eq. 
(p0|) . where the most general bi-local quantum map 
£[p] = £a® £b[p] operating upon the generic input two- 
qubit state p, is represented by a Kraus decomposition: 

ps =£ A ®E B \p] = Y J ^K{A tl ®B v ) p(Al®Bl), (51) 



where now the equality symbol can be used since we as- 
sume that both single-qubit maps £a and £3 are trace- 
preserving, 



3 3 



(52) 



i/=0 



but not necessarily unital: £f[I] ^ I, F £ {A, B} |3q . 
Under the action of £ , the initial state of each qubit trav- 
elling in path A or path B is transformed into either the 
output state 



p£ = Tr\ B (p £ ) 



(53) 



or 



p§ = Tr\ A (pe) = £ KB uP B Bl (54) 

respectively, where p A = Tr| b(p), and p B = Tr\ a{p)- 
Without loss of generality, we assume that the two qubits 
are initially prepared in the singlet state: p — p s . Then 

Eqs. §MB reduce t0 Ps = E«^A F G {A,B}. 
From the previous analysis [see Eqs. (|30H32p ] we know 
that to each bi-local quantum map £ a £3> £ b can be asso- 
ciated a pair of classical Mueller matrices Ma and Mb 
such that 



(pf) 



(55) 



The real-valued Mueller matrices Ma and M B associated 
via Eq. (fT2"|) to Ma and Mb, respectively, can be written 
as 



Ma = 



1 T 
a A 



M B = 



1 T 
b B 



where Eq. (fTS"]) with d 
been used, and 



A 



= ds and Mq 



(56) 



1 has 



(57) 



are the polarizance vectors of Ma and Mg , respectively. 
We remember that the condition = d# = is a 
consequence of the fact that both maps £ a and £b are 
trace-preserving, while the conditions a ^ and b ^ 
reflect the non- unital nature of £a and £b- With this 
notation we can rewrite Eqs. (|53H54|) as 





ai 




bi 


a = 


a 2 


b = 


b 2 




«3 




63 



1 3 

2 a A^A" 



(58) 
(59) 



9 



where we have defined clq = 1 = &o- Moreover, the output 
two-qubit density matrix p£ = £{p s ] can be decomposed 
into a real and an imaginary part as ps = 
where 



l Ps > 



< 


P+ 


7+ 


6+ 




at 


£_ 


7- 


7+ 


5- 


a- 


0- 


<h 


7- 


P- 


a_ 



(60) 



and 



with 







r -e+ - 


-»?+ - 


-T+ " 




1 


e+ 


-T- - 


-77_ 


P l t- 


" 4 




- 








- t+ V- 


£- 







(1 + 


03) ±[63(1 


+ a 3 ) 


— C33 


a± = 


(1- 


-as) ±[63(1 


- a 3 ) 


+ C33 



(61) 



(62) 



and 



6± 



P± = h ± (0361 - C 3 i), 
7± = 01 ± (aife 3 - C13), 
ai&i - C11 =F (0262 - C22). 



(63) 



and 



£± = 6 2 ± («3^2 - C32), 
ry± = a 2 ± (a 2 6 3 - C 23 ), 
t± = a 2 &i - C21 ± (ai&2 - C12), (64) 

where Cy = {AB T ) l3l i,j E {1,2,3}. 

At this point, our goal is to determine the two vectors 
a, b and the two 3x3 matrices A, B such that p £ m = 
and 



Re 

Pe = Pmems 



g(p)/2 p/2 

l-g(p) 



p/2 g(p)/2 



(65) 



where 



g{p) 



2/3, 



P, 



< p < 2/3, 
2/3 <p < 1. 



To this end, first we calculate a and b by imposing: 



Pi 



P B S 



A 

Pmems 



B 

Pmems 



1 - g(p)/2 
g( P )/2 

fl(p)/2 
l- fl (p)/2 



(66) 

(67) 
(68) 



respectively. Note that only fulfilling Eqs. (|67H68p , to- 
gether with pg c — pmems and p\ m = 0, will ensure the 
achievement of true MEMS states. It is surprising that 



in the current literature the importance of this point 
is neglected. Thus, by solving Eqs. (|67H68[) we obtain 
fli = a-2 = 0, 03 = 1 — g(p), and b = —a, where Eqs. 
(|58H59[) have been used. Then, after a little of algebra, 
it is not difficult to find that a possible bi-local map 
£ = £ a ® £b that generates a solution ps for the equa- 
tion ps — pmems, can be expressed as in Eqs. (|55H56j) in 
terms of the two real- valued Mueller matrices 



M A 



M R = 



1 



^/p 

l-g[p) g(p) 



1 

-y/p 

vp 

g(p)-l -g(p) 



(69) 



It is easy to check that both Ma and Mb are physically 
admissible Mueller matrices since the associated matri- 
ces Ha and Hb have the same spectrum made of non- 
negative eigenvalues {A M } = {Ao, Ai, A2, A3}. In partic- 
ular: 

{A M } = {0, 1-p, 0, 1+p}, for 2/3<p< 1.(70) 
and 



f 1 5-71 + 36^ 5 + 7T+36^ | 
{A„}-|0,-, - , - j, (71) 

for < p < 2/3. It is also easy to see that the map £ 
can be decomposed as in Eq. (|51j) in a Kraus sum with 
A = A 2 = 0, 



Ax-JX 



yH^P 





^3VA 3 



and B Q = B 2 = 0, 





S1VA1 



V^p 



B3 V A3 



1 

Vp 



-vp 

1 



(72) 



(73) 



for 2/3 < p < 1. Analogously, for < p < 2/3 we have 
^0 = 0, 



A2VX2 = 
and B = 0, 



A iy /X 



-<t>- 

v+ 



1/V3 




Bo v A2 = 



SiVAr = 



v+ 

0_ 



, A 3 VX 3 = 




1/V3 



V- 



A3 = 







(74) 



,(75) 



(76) 



,(77) 
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where 



/l / 1 + 6 P 
2 V VI + 36p 



(78) 
(79) 




1 ■ 



1 - 9p 



V1 + 36 P/ 

Note that these coefficients satisfy the following relations: 



3 



fi=0 



1 



(80) 



A straightforward calculation shows that the single-qubit 
maps £a and £ b are trace-preserving but not unital, since 



In this subsection, we show that both single-qubit 
maps £ a and £b can be physically implemented as lin- 
ear optical networks @ where polarization and spatial 
modes of photons are suitably coupled, without acting 
upon the pump beam. The basic building blocks of 
such networks are polarizing beam splitters (PBSs), half- 
waveplates (HWPs), and mirrors. Let \i,N) be a single- 
photon basis, where the indices i and N label polariza- 
tion and spatial modes of the electromagnetic field, re- 



spectively. We can also write \i,N) 



it 

iN 



|0) in terms of 



the annihilation operators &iN and the vacuum state |0). 
A polarizing beam splitter distributes horizontal (i = H) 
and vertical (i — V) polarization modes over two distinct 
spatial modes, say N — n and N — m, as follows: 



\H,n) in -> | if, n) ou t and |F,n) in -> \V,m) c 
\H,m)i a -> |if,m) ou t and \V,m} in -> \V,n) c 



(84) 



as illustrated in Fig. 3. A half-waveplate does not cou- 



and 



fj,=0 



v=0 



- g(p) o 
o g(p) 



g{p) 
o 2 



(82) 



(83) 



At this point our task has been fully accomplished. How- 
ever, before concluding this subsection, we want to point 
out that both maps £a and £b must depend on the same 
parameter p in order to generate proper MEMS states. 
This means that either a classical communication must 
be established between Ta and Tb in order to fix the 
same value of p for both devices, or a classical signal en- 
coding the information about the value of p must be sent 
towards both Ta and Tb . 



1. Physical implementation 

Now we furnish a straightforward physical implemen- 
tation for the quantum maps presented above. Up to 
now, several linear optical schemes generating MEMS 
states were proposed and experimentally tested. Kwiat 
and coworkers [38J] were the first to achieve MEMS using 
photon pairs from spontaneous parametric down conver- 
sion. Basically, they induced decoherence in SPDC pairs 
initially prepared in a pure entangled state by coupling 
polarization and frequency degrees of freedom of the pho- 
tons. At the same time, a somewhat different scheme was 
used by De Martini and coworkers [4l[ who instead used 
the spatial degrees of freedom of SPDC photons to in- 
duce decoherence. In such a scheme the use of spatial 
degrees of freedom of photons required the manipulation 
of not only the emitted SPDC photons, but also of the 
pump beam. 



I i,m) 
!■" 'out 



\i,n) 



-*■ \i,n. 



out 



J,m) 



FIG. 3: The polarizing beam splitter couples horizontal and 
vertical polarization modes (i, j € {H, V}), with two distinct 
spatial modes N = n and N — m of the electromagnetic field. 

pie polarization and spatial modes of the electromagnetic 
field and can be represented by a 2 x 2 Jones matrix 
Thwp(0) as 



Ti 



HWP 



(0) 



cos 26 -sin 29 
sin 26* cos 26 



(85) 



where 9 is the angle the optic axis makes with the hori- 
zontal polarization. Two half-waveplates in series con- 
stitute a polarization rotator represented by Tr(0) = 
Thwp(Qq+9/2)Thwp(Qq): where 8 is an arbitrary angle 
and 



Tr(9) 



cos 9 
sin 6 



— sin# 
cos 9 



(86) 



By combining these basic elements, composite devices 
may be built. Figures 4 (a-b) show the structure of a 
horizontal (a), and vertical (b) variable beam splitter, 
denoted HVBS and VVBS, respectively. HVBS performs 
the following transformation 



\H, n) 



cos 9\ H, n) 



\H, m) 



(87) 
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while VVBS makes 



\V,m) 



cos 9\V,n) 



■ singly, m) 



(88) 



At this point we have all the ingredients necessary to built 
the optical linear networks corresponding to our maps. 
We begin by illustrating in detail the optical network 
implementing £a (for 2/3 < p < 1), which is shown in 
Fig. 5. Let |^o) = a\H)+b\V) be the input single-photon 
state entering the network. If we define the VVBS angle 



9 p = arccos y/p, 



(89) 



then it is easy to obtain after a straightforward calcula- 
tion: 



\a{) = V%>i|$>) =b^T^p\H), 



\a{) = \/%A 3 \ipo) = a\H) + b^/p\V). 



(90) 



(91) 



Since detector Da does not distinguish spatial mode 1 
from spatial mode 2, the two states \a{) and \a$), sum 
incoherently and the single-photon output density matrix 
can be written as ps A = |a{)(a{| + |a|)(a!3|, where 



P£a 



\b\ 2 {l-p) ab*^p 



a*by/p 



p\b\< 



(92) 



Of course, if we write the input density matrix as po 
\i J o)(' l Po\: it is easy to see that 



P£a =J2 X ^P° A t- 

fj,=0 



(93) 




~7\ PBS \\ HWP 
VVBS I Mirror 



FIG. 5: Linear optical network implementing £a (for 2/3 < 
p < 1), for MEMS I generation. 



where Eqs. (J72J) have been used. Equation ([93]) . together 
with Eq. (|53|) . proves the equivalence between the quan- 
tum map £a and the linear optical setup shown in Fig. 
5. Note that the Mach-Zehnder interferometers present 
in Figs. 5 and 6 are balanced, that is their arms have the 
same optical length. In a similar manner, we can phys- 
ically implement £b (for 2/3 < p < 1), in the optical 
network shown in Fig. 6, where we have defined 



= V^iB^o) =b^T^p~\V), 



(94) 



(a) 



s'mO \H,m) 



HVBS= |tf."} in ^^~ 



■ cos6»#,. 



2 



(b) 



sin#|V,m 



VVBS : 



+^ 

a + 2 




cos e\V,n) 



PBS HWP 



FIG. 4: The variable beam splitters HVBS and VVBS. 



|/? 3 7 ) = v^Wo) = -bVP~\H) +a\V), 



(95) 



and, again, P£b = |#></3f| + |/3f)</3fl- 

The optical networks necessary to realize quantum 
maps generating MEMS II states are a bit more com- 
plicated. In order to illustrate them we need to define 
the following two angles 9 1 / 3 and 9^, that determine the 
transmission amplitudes of two VVBSs used in the MEMS 
II networks: 



(96) 
(97) 



In addition, a third angle 9^ determining the transmis- 
sion amplitudes of a HVBS, must be introduced: 




arccos < 



(98) 



Then, the map £a (for < p < 2/3), is realized by the 
optical network shown in Fig. 7, where we have defined 



K J ) = V%4#o) = -a^ \H) +bi>+\V), (99) 
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▲ 



V 



]45° 



^| PBS [] HWP 
VVBS 1 Mirror 



FIG. 6: Linear optical network implementing £b (for 2/3 < 
P < 1), for MEMS I generation. 




\/\ PBS \] HWP 



VVBS 



HVBS | Mirror 



\*") = V^^slV'o) = a0 + |JT>+ 6^-1*0. (100) 



K 7 ) = v ^A 1 |^o)-^|^>. (101) 

In this case, incoherent detection produces the output 
mixed state p £A = \a^)(a^\ + |a^ 7 )(a| / | + |a( / )(a{ / |. 
Finally, the map £b (for < p < 2/3), is realized by the 
optical network shown in Fig. 8, where we have defined 

\ti I ) = ^B 2 \i/j )=biP + \H)+a<j>-\V), (102) 



|/?3 77 } = y/^B 3 \^o) = |#) + (103) 



l/?i /7 > - v/AT^il^o) - ^l^>- (104) 
As before, now we have ps B = + + 



VI. SUMMARY AND CONCLUSIONS 

Classical polarization optics and quantum mechanics 
of two-level systems are two different branches of physics 
that share the same mathematical machinery. In this 
paper we have described the analogies and connections 
between these two subjects. In particular, after a review 
of the matrix formalism of classical polarization optics, 
we established the exact relation between one- and two- 
qubit quantum maps and classical description of linear 
optical processes. Finally, we successfully applied the 
formalism just developed, to two cases of practical utility. 



FIG. 7: Linear optical network implementing £a (for < p < 
2/3), for MEMS II generation. Each of the two Mach-Zehnder 
interferometers constituting the network are balanced. 




\Yo) + 



^ PBS [] HWP 
VVBS 1 HVBS I Mirror 



FIG. 8: Linear optical network implementing £b (for < p < 
2/3), for MEMS II generation. Each of the two Mach-Zehnder 
interferometers constituting the network are balanced. 



We believe that the present paper will be useful to 
both the classical and the quantum optics community 
since it enlightens and puts on a rigorous basis, the so- 
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widely used relations between classical polarization op- Acknowledgments 
tics and quantum mechanics of qubits. A particularly 
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